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ABSTRACT 

The Advanced Amphibious Assault Vehicle (AAAV) is a high water speed 
amphibious armored personnel carrier that will replace the current family of Marine 
Corps amphibious assault vehicles. The AAAV is currently in Phase I of the DOD 
Acquisition Process. During this phase-extensive development, testing is being 
conducted and prototypes are being constructed. Ongoing tests of the current electrical 
system architecture are revealing problems and issues that need to be addressed. Present 
testing is also revealing the need for a troubleshooting tool that can be used to simulate 
and test proposed solutions. 

A distributive computing effort is being conducted with Purdue University in 
order to provide the Marine Corps a tool where by it can test and evaluate the entire 
AAAV electrical system architecture. This document provides a general system 
description of the AAAV, an electrical system architecture overview, and a candidate 
electrical system description. Testing, modeling, and computer simulation efforts applied 
to the NBC (Nuclear, Biological, and Chemical) fan/filter motor drive and battery banks 
are described and the resulting data analyzed. Current research is investigating an issue 
regarding energy being forced back into the high power distribution bus due to 
regeneration from sudden impacts on the vehicle turret. An ultra capacitor bank may be 
incorporated in the AAAV electrical system in an effort to mitigate the effect of this 
regenerative energy. Therefore, testing, modeling, and computer simulation of an ultra 
capacitor is also reported. 
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I. 



INTRODUCTION 



A. AAAV Program 

The AAAV (Advanced Amphibious Assault Vehicle) will provide the Marine 
Corps a weapons system fully capable of implementing ship-to-objective maneuver as an 
integral part of the amphibious triad (AAAV, V-22, LCAC) to execute the concepts of 
Operational Maneuver From The Sea (OMFTS) and Ship to Objective Maneuver 
(STOM). Battlespace dominance by Marine forces will be significantly enhanced as a 
result of the AAAV’s high water speed and superior land mobility which have 
historically limited the rapid maneuver of armored combat vehicles. The AAAV will 
travel as fast as 29 miles-per-hour (mph) on water as opposed to the 10 mph the current 
AAV can travel. The AAAV is designed to allow immediate, high-speed maneuver of 
Marine infantry units as they emerge from attack positions aboard ships located beyond 
the visual horizon -- 25 miles and beyond. Projection of these forces will be conducted 
as a single, seamless stroke that capitalizes on the intervening sea and land terrain to 
achieve surprise and rapidly exploit weak points in enemy littoral defenses. In short, the 
AAAV is faster, can travel greater distances before refueling, is navigationally smarter 
and more maneuverable than its predecessor. 

The AAAV is currently in Phase I (Program Definition and Risk Reduction) of 
the Department of Defense (DoD) Acquisition Process. Phase II (Engineering and 
Manufacturing Development) is scheduled to begin Fiscal Year 2001 (FY01). The final 
phase, Phase III (Production, Fielding/Deployment and Organizational Support), is 
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scheduled to begin FY06. The total cost of the program is expected to be $7.5B with 
each unit costing S6.6M. The Marine Corps is scheduled to receive 1013 units. [Ref.l] 

B. Research Focus 

A distributive computing effort is being conducted with Purdue University in 
order to provide the Marine Corps a tool where by it can test and evaluate the entire 
AAAV electrical system architecture. Reporting in this thesis includes a general system 
description of the AAAV, an electrical system architecture overview, and a candidate 
electrical system description. Efforts are then applied to documenting the testing, 
modeling, and computer simulation of the NBC (Nuclear, Biological, and Chemical) 
fan/filter motor drive and battery banks. Furthermore, research will also address the 
requirement to capture energy being forced back into the high power distribution unit due 
to regeneration from sudden impacts on the vehicle turret. An ultra capacitor bank may 
be incorporated in the AAAV electrical system in an effort to capture regenerative 
energy. Therefore, testing, modeling, and computer simulation of an ultra capacitor is 
also reported. 

C. Chapter Overview 

Chapter I of this thesis provides an introduction to the AAAV program. It also 
identifies the research focus and provides a chapter overview. Chapter II contains a 
comprehensive system description of the AAAV. It details a general system description 
and introduces the four modes of operation. In addition the discussion within Chapter II 
addresses the electrical system architecture, beginning with a topology overview and 
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ending with a primary components description. Chapter III contains an introduction to 
distributed computing and background information on the programming language ACSL 
(Advanced Continuous Simulation Language), which will be used in this thesis to 
develop an NBC fan drive model. Further, the Automated State Model Generation 
Algorithm used by P. C. Krause & Associates (PCK&A), and Purdue University, to 
develop their NBC fan drive model is described and referenced in Chapter III. Chapter 
IV includes a description of the selected candidate electrical system architecture. In 
addition it provides background on the major components of the selected architecture, 
specifically the NBC fan, ultra capacitor, and batteries. Chapter V provides discussion 
regarding the modeling efforts applied to the NBC fan drive, ultra capacitor, and battery 
unit. It also provides results and conclusions for each model. Chapter VI is the 
conclusion. It summarizes the research while documenting results of the work. It also 
identifies Department of Defense relevance of the work and ideas for future work. 
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II. AAAV SYSTEM DESCRIPTION 



A. Introduction 

This chapter provides a comprehensive system description of the AAAV. It 
details a general system description and introduces the four modes of operation. In 
addition the chapter addresses the electrical system architecture, beginning with a 
topology overview and ending with a brief description of the primary components. 

B. System Description 

The Advanced Amphibious Assault Vehicle (AAAV) is a high water speed 
amphibious armored personnel carrier that will replace the current family of Marine 
Corps amphibious assault vehicles, the AAV7A1 series. Operationally configured, the 
AAAV will weigh about 37 tons, be able to carry 17 combat-equipped Marines and a 
crew of three over 3 foot high waves in excess of 20 knots, and travel over land at 45 
mph. Armed with a 7.62 mm machine gun and a 30 mm cannon, the AAAV will use 
Global Positioning System, forward-looking infrared, and a night vision system for 
navigation, targeting, and intelligence gathering. [Ref. 1] 

The AAAV force is designed to provide a tactical assault and sustainment 
capability for Marines of the Marine Air-Ground Task Force. From amphibious ships 
standing well offshore, even over the horizon from the objective, the AAAV rapidly 
transports the landing force over the beachhead to an objective ashore, using maneuver 
and speed, plus on-board firepower to achieve superiority over enemy forces. Once 
ashore, the AAAV will serve as an armored personnel carrier, providing transportation. 
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protection, direct fire support and command, control and communications (C3I) for the 
Marines. 

High water speed requires minimal hydrodynamic drag; to accomplish this, the 
AAAV will retract its tracks and suspension system and deploy cover plates on its 
underside to present a smooth surface. Thus, while in the water and traveling at high 
speed, the AAAV will be a planing hull craft. When the AAAV approaches 15 to 20 feet 
of water depth at the shore, it will reduce speed to come off plane, retract the cover 
plates, lower the suspension and tracks, and travel the remaining distance to the shore at a 
speed of approximately 8 to 10 knots. 

Once ashore, the Marine Corps will use the AAAV (P) as a tracked armored 
personnel carrier to accomplish dominant maneuver and precision engagement', in that 
mode, the 1 7 infantrymen carried inside will typically dismount to fight. The crew will 
use the AAAV (P)’s primary and secondary weapon (predicted to be a 25-35 mm cannon 
and 7.62 mm coaxial mounted machine-gun, respectively, with fully stabilized turret and 
forward looking infrared optics) to support the infantry and armor combat elements. 

The AAAV’s land mobility characteristics must be comparable to the Marine 
Corps’s Ml A1 Abrams main battle tank. This requires a top speed of approximately 45 
mph, the capability to traverse the same terrain as the tank during cross-country 
operations, and the capability to cross the same obstacles and terrain features (for 
example, trenches, hills, walls, and soft soils) as the tank. 
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Figure 2-1, Advanced Amphibious Assault Vehicle. [Ref. 1] 

1. System Operating Modes 

a. High Speed Water Mode 

An exploded view of the AAAV in high speed water mode is 
provided in Figure 2-2. The AAAV will operate in this mode 
during waterborne operations and will possess the following 
attributes. 

• Top Speed: 20-25 knots in 3 feet significant wave height 

• Maintains a heading in 1 0 feet significant wave height 

• Self-righting ( 1 00 degrees) 

• High water speed turning diameter <170 meters 

• Stops forward motion within 75 meters 

• Water range: 65 nautical miles 
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Extended Bow FI 



Deployed Chine Fla 



Retracted Suspension & Track 



Power Transfer Module 



Deployed Transom Flap 



Figure 2-2, AAAV High Speed Water Mode. [Ref. 1] 

b. Transition Mode 

A view of the vehicle in transition mode is documented in Figure 
2-3. The AAAV will operate in this mode while entering and 
exiting the water (from ship to land). Additional details include 

• Used to enter and exit the water (from ship to land) 

• Slow waterspeed operations 

- Top speed: 9 knots (4.5 knots reverse) 

Surf transit: 8-foot plunging surf 

- Survivability: 360° rollover, 180° pitch-poling 

- Fully submersible for 5 seconds 

• Also used for riverine and coastal operations 
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Forward & Reverse 
Steering 



Power Transfer Module 



Figure 2-3, AAAV Transition Mode. [Ref. 1] 

c. Land Mode 

The vehicle in land mode is pictured in Figure 2-4. The AAAV 
will operate in this mode on land and possess the following 
attributes. 

• Cross country mobility: > = Ml A1 Abrams tank 

• Top speed: 45 mph (0-20 mph in 7 seconds) 

• Range: 400 miles (land only) or 250 miles after 1-hour (25 
nautical miles) high water speed transit 

• Obstacle crossing 

8-foot trench span, 3-foot vertical wall 
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Figure 2-4, A AAV Land Mode. [Ref. 1] 

d. Silent Watch Mode 

A number of the key system components associated with the silent 
watch mode are illustrated in Figure 2-5. The AAAV will operate 
in this mode while stationary on land. In this mode, the following 
attributes hold. 

• Used for stationary vehicle to reduce acoustic and thermal 
signature 

• Fully combat ready 

• One or more crewman on watch 
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All com/nav/fire control/vision subsystems fully 
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Figure 2-5, AAAV Silent Watch Mode. [Ref. 1] 



C. Electrical System Architecture 
1. Topology Overview 

A block diagram of the electrical system is shown in Figure 2-6. This architecture 
schematic is of a level of detail greater than will be initially modeled by P.C. Krause & 
Associates or addressed in this thesis. The candidate electrical system for modeling and 
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simulation can be found in Chapter IV, Figure 4-1 . While a detailed description of the 
entire electrical architecture is outside the scope of this thesis, the primary components 
are identified below. 




Figure 2-6, Electrical System Architecture. [Ref. 2] 
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2 . 



Primary’ Components 



a. Batteries 

- Eight (8), 12 volt, manufactured by Optima 

• Three series pairs dedicated for starting. 

• One series pair dedicated to the Vetronics Bus during starting 

b. Generators 

• Main generator - 800 Amp, six-phase wound-rotor 
synchronous generator connected to two six-pulse diode bridge 
rectifiers whose dc outputs are connected in parallel. The field 
windings are connected to an exciter that is a three-phase 
machine whose armature voltages are rectified using a rotating 
six-pulse diode bridge rectifier. 

• Auxiliary Generator - 330 Amp, connected to the 
Vetronics Bus 

c. Ultra Capacitor 

• 100 Volt, IF, 5000J capacitor manufactured by Pinnacle 
Research Institute. Not shown in Figure 2-6. See Figure 4-1. 

d. Hull Power Distribution Unit (HPDU) 

• Provides main power distribution to the Vetronics bus and high 
power bus. 
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e. Remote Acquisition and Control Modules (RACMs) 

• Five (5-HRACM) hull and one (1-TRACM) turret 

• Provide Secondary power distribution, Sensor Data Collection, 
Sensor Power, and Communications through Utility Data Bus 

f. Sealed Slip Rings 

• Provides power and data signal path to electrical systems 
located in the turret. 

g. NBC Fan Drive 

• 3.5 Hp, Wye wound, permanent magnet, brushless DC motor 
manufactured by Invincible Airflow Systems. Powered using a 
120-degree voltage source inverter with sensorless control. 
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III. COMPUTER SIMULATION 



A. Introduction 

Chapter III contains an introduction to the theory behind the application of 
distributed computing to power system simulations. It also provides an overview of the 
programming language ACSL (Advanced Continuous Simulation Language), which will 
be used in this thesis to develop an NBC fan drive model. Further, it will introduce the 
Automated State Model Generation Algorithm used by PCK&A, and Purdue University, 
to develop their NBC fan drive model. 

B. Distributed Computing 

The consulting firm, P.C. Krause and Associates (PCK&A), has developed a 
technique that will aid in the simulation of complex systems in general and power 
electronic based systems in particular. This technique has been termed distributed 
computing. Distributed computing is a heterogeneous simulation environment that 
makes it possible to divide the simulation of power electronic based systems into any 
number of synchronized personal computers (PC). Moreover, subsystems can be 
implemented using different languages (i.e., Simulink, Saber, and ACSL) making it 
possible to have the advantage of a very large computer by paralleling inexpensive PC’s. 

In addition to the gain in computational speed, another advantage of this 
simulation environment is the ability to decouple the user interface and the simulation 
that is controlled by the interface. In particular, the user interface can be executed on any 
computer that is connected to the Internet allowing users at remote locations to “exercise” 
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the simulation and retrieve data. Also vendors who develop proprietary simulations can 
link the model as executable and not divulge the model inner workings. 

C. ACSL Background 

The acronym ACSL (pronounced “axel”) stands for the Advanced Continuous 
Simulation Language. It has been developed expressly for the purpose of modeling 
dynamic systems described by time dependent, nonlinear differential equations, and/or 
transfer functions [Ref. 3]. Simulation of physical systems is a standard analysis tool 
used in the evaluation of hardware design prior to and in concert with actual construction. 
An accurate simulation model may provide invaluable insight into design flaws or system 
limitations before extensive capital investment is made on hardware. For instance in 
analyzing the AAAV electrical power system, the engineers of General Dynamics 
uncovered an instability in the electrical system during low power operation when the 
batteries are disconnected. Discovering the instability early will allow the General 
Dynamics team to troubleshoot the instability. Knowing this up front, the proposed 
design could then be modified without costly retrofits or adhoc solutions. 

For the power engineer, ACSL facilitates the analysis and control of complex 
electromechanical networks. Variable speed electric drives are popularly used in • 
pumping applications (fans, blowers, and compressors), large machinery actuation 
(elevators, hoists, cranes, conveyors, and crushers), remote servoing (robotics and 
precision machine tools), electric vehicle propulsion (cars, busses, trains, and ships), and 
as an alternative to hydraulic and pneumatic actuation systems. In addition to the 
electromechanical device, all of these systems are characterized by nonlinearities which 
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make application of linear system theory difficult (if not impossible). After designing 
comprehensive electric drives containing power electronic converters, electric machines, 
multi-loop analog and digital controls, and various protection and sensing devices, it is 
advantageous for the engineer to develop computer-based models which predict as 
accurately as possible the operational modes of the interconnected system. As the 
hardware becomes available, the models are validated and refined as necessary. Thus, 
modeling and simulation are integral parts of the design and analysis process. 

The ACSL program is intended to provide a simple method of representing these 
mathematical models on a digital computer. Working from an equation description of the 
problem (with logical conditions conveniently handled), the user writes ACSL statements 
to describe the system under investigation. Run-time commands exercise the model; 
these statements can be submitted for solution in command files or entered interactively. 
Simulation studies may be conducted, variables plotted, and the system re-evaluated for 
changes in parameters and operating conditions. 

D. Automated State Model Generation Algorithm 

P. C. Krause and Associates (PCK&A) are using a new approach of automatically 
generating state space models in their simulation, and study, of the candidate electrical 
system for the AAAV. This new approach is termed the Automated State Model - 
Generation Algorithm [Ref. 4], As the software is proprietary and not widely available, 
this method will not be used in the simulation and study conducted by this thesis. This 
background is included to ensure a complete and thorough understanding of the methods 
employed by PCK&A in their development and study of the candidate electrical system. 
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background is included to ensure a complete and thorough understanding of the methods 
employed by PCK&A in their development and study of the candidate electrical system. 

In the Automated State Model Generation Algorithm the candidate electrical 
system is described by the standard node incidence matrix and elementary branch data 
(e.g. resistances inductances, and back emf s) [Ref. 4], The standard node incidence 
matrix is established by numbering the nodes and branches of the circuit. For example in 
Figure 3-1 below, the node incidence matrix would be of the form shown in Figure 3-2. 
The nodes are labeled vertically on the right while the branch elements are specified 
horizontally across the top. 




Figure 3-1, Example Circuit. [Ref. 4] 
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Figure 3-2, Node Incidence Matrix. [Ref. 4] 



The elementary branch that PCK & A employs to model portions of systems ranging 
from thyristor switches to machine phases is illustrated in Figure 3-3. If the switch is closed 
(state=l), the branch is active. If the switch is open (state=0), the branch is inactive. 



State = 0 e { Ci 




Figure 3-3, Branch Model. [Ref. 4] 



From this input data the Automated State Model Generation Algorithm will automatically 
establish the composite system state equations, thus relieving the analyst from establishing the 
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differential equations by more tedious conventional means [Ref. 4]. Since the state equations are 
now available, linearization and eigenvalue analysis can be performed at the system level. 

Further, the state equations can then be solved numerically using one of a number of well- 
established numerical algorithms (i.e., ones available in ACSL: Runge-Kutta, Gears-Stiff, etc.). 
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IV. ELECTRICAL SYSTEM COMPONENT BACKGROUND 



A. Introduction 

This Chapter includes a description of the selected candidate electrical 
system architecture. In addition it provides background on the major components of the 
selected architecture, specifically the NBC fan drive, ultra capacitor, and batteries. 

B. Candidate Electrical System Description 

Figure 2-6, illustrates a general electrical system architecture for the 
AAAV. From this system architecture PCK&A selected a candidate system, Figure 4-1, 
for initial modeling purposes and to demonstrate the principle of distributed computing. 
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Main Auxiliary 

Generator Generator 




Main APU 

Engine Engine 

Starter Starter 



Figure 4-1, AAAV Candidate Electrical System Description. [Ref. 2] 

C. NBC Fan Drive Background 

1. Introduction 

The brushless dc motor is becoming widely used as a small-horsepower control 
motor. The brushless dc motor, however, is really not a dc motor at all; in fact it is a 
permanent-magnet synchronous machine. Nevertheless, when it is supplied from a 
source with frequency equal to the electrical angular speed of its rotor, it becomes a 
brushless dc motor. The nomenclature is adopted not because it looks like a dc motor, but 
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because its torque speed characteristics can be made to resemble those of a dc shunt 
motor with a constant field current. If the permanent-magnet synchronous machine is to 
be supplied from a source of frequency that always corresponds to the electrical angular 
speed of its rotor, it is necessary that the rotor position be identified. This can be done in 
several ways, such as with a light beam and a toothed gear, or indirectly by using the 
back emf voltages found by sensing the open-circuit phase voltages, or by using Hall- 
effect sensors which magnetically sense the position of the rotor poles. The stator voltage 
source and rotor electrical angular speed are self-synchronized using a dc-to-ac inverter, 
wherein the transistors are switched on and off based on rotor position. There are 
numerous types of inverters. This thesis will explore the 180° and the 120° voltage source 
inverters. 

In this section a description of the permanent-magnet synchronous machine is set 
forth. Next, the 180° inverter operation is explained. The operation of the 180° inverter 
is then used as a basis for the explanation of the operation of the 120° inverter. This 
explanation is completed by defining a sensorless control method for the 120° inverter. 

2. Description of the Brushless DC Machine 

The term "brushless dc" is somewhat of a misnomer since the brushless dc 
machine is actually a permanent-magnet synchronous motor such as the one depicted in 
Figure 4-2. In fact, the only difference between the common notion of a permanent- 
magnet synchronous machine and the brushless dc machine is in the form of the applied 
voltages. 
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In this chapter a background on the brushless dc machine is provided beginning 
with a description of the permanent-magnet synchronous machine in machine variables. 
The machine equations are then transformed into the rotor reference frame. Finally, in 
the context of the rotor reference frame, the form of the applied voltages is specified. 
Additional background is available in references [5 ] and [6], 



bs axis 




as axis 



Figure 4-2, Three-Phase Permanent-Magnet Synchronous Machine. 



a. Machine Equations 

The purpose of this section is to set forth the machine equations for the 
permanent-magnet synchronous machine. The equations presented here are subject to the 
following assumptions: (1) the machine is non-salient, (2) the back emf is sinusoidal, and 
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(3) the machine is wye connected. The stator phase voltage equations describing the 
permanent-magnet synchronous machine are given in vector form by (4.1) 



abcs 



= r a 



s abcs 



+ pA 



abcs 



(4.1) 



where p is the Heaviside operator d/dt. The flux linkage equations are similarly given by 
(4.2) 
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L -M 
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(4.2) 



where X m is a machine parameter dependent on the strength of the permanent magnet. 

Since the machine is wye connected ( the three stator currents must sum to 
zero), (4.2) may be rewritten as (4.3) where L s is defined by (4.4). This representation is 
convenient since -the three stator equations become uncoupled. 



^ abcs ~ L s i a bcs + 



sin ( 0 r ) 

• (a 2 *0 

l 3 



J 



*T' + f, 



(4.3) 
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L S =L + M 



(4.4) 



The electromagnetic torque ( T e ) is established by expressing the co- 
energy, taking the derivative of the co-energy with respect to d r , and multiplying the 
result by the number of pole pairs [Ref. 5], 



T e = 
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Jj 


m 

V 
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(4.5) 



b. Machine Equations in the Rotor Reference Frame 

Considerable insight and simplification may be gained by making a 
change of variables which transforms all machine variables (abcs) to rotor reference 
frame variables (qdOs). The transformation from abcs variables to qdOs variables is given 

by (4.6) and (4.7) where / may represent voltage, current, or flux linkage. 

/*„ = K; /*, (4.6) 
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The angle 6 r is the rotor electrical angle that necessarily varies at the rotor 



electrical angular velocity. The inverse transformation is 



ter - 



COS 



M sin(0 f ) 1 



cos 
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° r 3 J S 1 I V r 3 

2k^ 



0 r + — 

V 3 j 



sin 



J 

2k ^ 



d r + — 

\ 3 j 



(4.8) 



Applying the transformation to the machine equations and expanding 
those equations yields 

V qs = i r s + P L s ) iqs + op LA + co r A'n, (4.9) 



v * =( r s + pLs)i r ds-0) r L s i r gs 



(4.10) 



T = 






(4.11) 



It can be seen from (4.1 1) that if torque is to be constant, then i' s must be 

constant. Hence from (4.9) and (4.10), steady-state operation requires that v' and 
are constant as well. The steady-state line-to-neutral voltages may be given by 
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(4.12) 



V s cos(<? r + a) 




(n 2 * ^ 


V' cos 
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l 3 J 
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V s cos 


6 r + — + a 



where a is some arbitrary phase angle. Note, to achieve brushless dc operation, the stator 
voltages must vary at the rotor electrical angular velocity. Then, for any rotor speed, 

(4.13) and (4.14) apply. If the amplitude V s is constant, then so are v r qs and v r ds . 



v g S = K cos (a) 



(4.13) 



v * = ~K sin(a) 



(4.14) 



For steady-state conditions, with a = 0 and co r L s « r s , (4.13), (4.14), 
(4.9)-(4.1 1) may be manipulated to yield (4.15) which is identical in form to the steady- 
state torque equation of an armature-controlled dc machine. 







(4.15) 
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Now that the basic features of brushless dc operation have been set forth, 
the problem of actually obtaining the voltages described by (4.12) is addressed. One 
device used to obtain these voltages in practice is known as an inverter. The inverter is 
the subject of the next section. 

3. Inverter Operation 

There are several methods of achieving stator voltages of the form given by 
(4.12), where the frequency of the voltages is variable. One method is to use a device 
known as an inverter. A typical three-phase inverter is depicted in Figure 4-3. Assuming 
ideal components, the inverter operates without loss and thus a conducting device is 
modeled as a short circuit and a blocking device is modeled as an open circuit. The 
tradeoff in using an inverter is that the voltages produced contain harmonics that increase 
loss in the stator circuits and can introduce torque pulsations. 

There are two types of inverters, the 120° and the 180°. In the 120° case, each 
transistor is on 120° out of a 360° cycle. Likewise, in the 180° case, each transistor is on 
1 80° out of 360°. Although the 120° inverter-driven brushless dc machine is more 
difficult to handle analytically, it has a practical advantage in that it can be operated 
without any rotor position sensors. In this case, position information is derived indirectly 
from measurement of the open-circuit phase voltages. 
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a. Basic Inverter Operation 



A schematic depiction of the brushless DC inverter system appears in Fig 
4-3. The inverter consists of six transistors, T1-T6, which are used to convert DC power 
to AC power and six free-wheeling diodes, D1-D6, which are necessary to protect the 
transistors from the inductive nature of the load. 




The first step toward the analysis of the inverter is to apply Kirchhoff s 
Voltage Law (KVL) starting at the negative dc rail (g), proceeding across the lower 
transistor and across the machine phase, and returning to the negative dc rail via the 
neutral point (n). Applying this path for each machine phase yields 
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(4.16) 






= v M + V, 



”g 



V bg 



= V bs + V - 



ng 



(4.17) 



eg 



= V + V 
cs r v ng 



(4.18) 



Kirchhoff s Current Law (KCL) dictates for a wye connection that 



L + hs + L = 0 (4-19) 

Equations (4.19) and (4.3) may be used to show that the sum of the stator flux linkages, 
namely A as + A bs + A cs , is zero. Since the sum of the stator currents is zero, and the sum 
of the stator flux linkages is zero, (4.1) shows that 



V +v, +v 

as bs 1 c 



= 0 



(4.20) 



The neutral to ground (where ground is the lower rail of the dc supply) voltage, v ng , is 
obtained by adding (4.16), (4.17), and (4.18), and then subtracting (4.20) which results in 



v = — (v + v. + v ) 

v ng ^ V a Z b g eg ) 



(4.21) 
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Substituting (4.21) into (4.16), (4.17), and (4.18) yields the following expressions for the 
phase voltages. 



2 1 1 

V as ~ n V °g ^ 3 Vcg 



(4.22) 
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(4.23) 
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b. Six-Step 180° Continuous-Current Inverter 

The logical signals for the switching (commutation) of the transistors are 
shown in Figure 4.4. Ti is turned on at 6 r = -90° and turned off at 6 r = +90° , at the 
same time T 4 is turned on. Note, in practice some delay must be present to avoid short- 
circuiting the dc input supply. The gating signals for the top transistors associated with 
the remaining two phases (T3 and T5) are each phase displaced from the Tj signal by 120° 
so that a balanced set of voltages is produced across the load. It is assumed that the turn- 
off time of the transistors is negligible whereupon the transistors become ideal switches. 
At the instant T 1 is turned off, and T 4 is gated on, the current T 1 was carrying is diverted 
to the diode D4. The diode continues to conduct until the current decreases to zero. Once 
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the current i as reverses direction it is carried by T4. This mode of inverter operation is 
generally referred to as the continuous-voltage or continuous-current mode. [Ref. 7] 



T1 
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Figure 4-4, Transistor Switching Logic for the 180° 
Continuous-Current Inverter. 
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The operation of the 180° inverter is set forth in Table 4.1. In the table 
electrical rotor position is divided into six 60° segments known as Basic Switching. 
Intervals (BSI). The second column specifies the center of each switching interval in 
terms of 0 r +<p. The angle <f> is known as the advance-firing angle, which allows the 
shifting of the transistor firing sequence relative to the position of the rotor. The third 
column of Table 4.1 defines the transistor firing sequence by indicating which transistors 
are gated on in a given BSI. The transistor firing pattern results in pole (abcg) and phase 
(abcs) voltages which are specified in Table 4.1 [Ref. 7]. Note that once the switching 
pattern is assigned, the phase voltages are completely determined. This fact makes 
programming the 180° inverter rather straightforward in ACSL. 
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Table 4.1, 180° Inverter Operation. 




Expressing the phase voltages given in Table 4.1 as a Fourier series in 
electrical rotor position, neglecting all harmonics, and incorporating the transformation to 
the rotor reference frame yields 



»; = ^cosW (4.25) 

n 

Vi — sin(^) (4.26) 

n 

Comparing (4.24) and (4.25) to (4.13) and (4.14), it is apparent that the 180° inverter a 
and (f) are equivalent. 



c. 120° Discontinuous-Current Inverter 

The logical signals for the switching (commutation) of the transistors are 
shown in Figure 4-5. Ti is turned on at 9 r = 0° and turned off at 6 r = +120° ; however, 
T 4 is not turned on until 60° later at 9 r = 180° Hence, the diode in parallel with T 4 begins 
conducting when Tj is turned off and continues to conduct until the current i as becomes 
zero. At this point the as winding of the machine becomes open- circuited. The 60° delay 
in turning T 4 on is done primarily for the purpose of lessening the possibility of a short 
circuit occurring through T 1 and T 4 and to provide an open-circuit period for indirectly 
sensing the rotor position. The voltage v ag in Figure 4-3 is equal to the source voltage 
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v dc if Tj is turned on or if the diode in parallel with T i is conducting. When T 4 is turned 



on or when the diode in parallel with T 4 is conducting, v ag is zero. 



Ti J 

T 2 | 

T 3 1 

T 4 1 

T 5 I ~ 

T 6 1 

0 1/2 x x 3/2 x 2 x 

Figure 4-5. Transistor Switching Logic for the 120° 
Discontinuous-Current Inverter. 

The operation of a 1 20° inverter is summarized in Table 4.2, which is 
much the same as Table 4.1 except for two important differences. The first is that only 
two transistors are gated on at any given instant of time. Secondly, question marks 
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appear periodically in the last three columns to indicate that the instantaneous current 
determines the voltages in the corresponding phase. [Ref. 7] 
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Table 4.2, 120° Inverter Operation. 

The dependence of the phase voltages on the current in the open-circuited 
phase is best illustrated by examining one specific BSI. In BSI II if i bs < 0, then current 
will flow through D3 and thus v bg = V dc . Conversely, if i b > 0 then current flows 

through D6 which causes v bg = 0. Finally, if i bs = 0, then v bs is the open-circuit voltage 
as derived from (4.1) and (4.3) and given by 



v 4* = cos 



e r 




(4.27) 



If (4.17) is substituted into (4.21), then the neutral to ground voltage is given by 
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(4.28) 



V = — (v + Vt + V ) 
ng 2 ' ag bs cg ' 



Upon using (4.27) for v bs in (4.28) and then using (4.28) to replace v ng in (4.17), the 
voltage across the lower transistor in leg B CU) is 



1 T/ 3,, 

''t g = 2 V * + 2 A ” a> r c os 






(4.29) 



Similarly, (4.28) can be used to find v as and v cs from (4.16) and (4.18). 
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This analysis can be repeated and adapted for each basic switching 
interval (BSI). Table 4.3 reflects the pole voltages, v abcg , when each phase current, i abcs , is 
exactly equal to zero. For each phase current, i abcs , greater than zero, each pole voltage, 
v abcg , is equal to V dc . For each phase current, i abcs , less than zero, each pole 
voltage, v abcg , is equal to zero. Table 4.4 similarly reflects the stator phase voltages, 
v abcs , when the phase currents, i abcs , are equal to zero. 
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Table 4.3, Pole Voltages for BSI I through VI with Non-Conducting 

Phase Open Circuited. 
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Table 4.4 Stator Phase Voltages for BSI I through VI with 
Non-Conducting Phase Open Circuited. 
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d. Sensorless Control 



Brushless dc motors are usually controlled using dc-ac inverters as 
previously discussed. The motor requires a rotor position sensor for starting and for 
providing the proper commutation sequence to turn on the power devices in the inverter. 
Hall effect sensors can be used; however, they must be specially mounted in the interior 
of the stator. Further these sensors are temperature sensitive, limiting the operation of the 
motor to below about 75 °C. Due to these limitations of the motor operation with Hall 
devices, sensorless operation of the PM brushless motor is receiving more attention. 

Several schemes for position sensorless operation of PM motors have been 
reported and reviewed. They are based on the following techniques. [Ref. 8] 

1. Position sensing using the back emf of the motor. 

2. Position information using the stator third harmonic component. 

3. Position information based on detection of the conducting interval of 
free-wheeling diodes connected in anti-parallel with the power 

devices. 

4. Position information using the motor parameters, and the measured 
terminal voltages and currents. 

Position sensing using the back emf of the motor is being utilized for 
control of the NBC fan motor. The sensed back emf can be used in several ways to 
obtain the proper switching sequence for the inverter power devices. Some or these 
methods are: 
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• The zero-crossing approach based on detecting the instant at which the 
back emf in the unexcited phase crosses zero or reaches a predetermined 
level. 

• Phase-lock-loop (PLL) technique to lock on to the back emf waveform in 
the unexcited phase winding during each 60 degree excitation interval. 

• Back emf integration approach, which provides a commutation pulse when 
the integrated rectified back emf waveform reaches a preset threshold 
value. 

The zero-crossing approach method is utilized in the NBC fan drive and therefore is 
discussed. In a 120 ° inverter-driven permanent-magnet brushless dc machine only two of 
the three phases are excited at any one time. The back emf voltage in the unexcited phase 
can be measured to establish a switching sequence for commutation of the power devices 
in the inverter. Based on the rotor position, the power devices are commutated 
sequentially every 60 degrees to continually synchronize the phase excitation with the 
rotating rotor magnetic field. 

The pole voltages v ag , v bg , v cg , and the neutral voltage v ng are measured. 
At the instants of the zero crossing of the back emf waveform, the pole voltage is equal to 

V 

the neutral voltage. When the back emf of phase A crosses zero, then v ag = v ng « . In 

order to use this zero-crossing point to derive the switching sequence, this point has to be 
phase shifted by 90 degrees [Ref. 8]. 
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D. Ultra Capacitor Background 

As shown in the candidate electrical system description, an ultra capacitor bank 
will be modeled and validated. This ultra capacitor may be incorporated in the AAAV 
electrical system in an effort to capture regenerative energy. Regenerative energy could 
be produced as a result of a sudden impact on the vehicle’s main gun turret. This energy 
could be forced back into the high power distribution bus thus potentially damaging the 
system. Theoretically, the ultra capacitor bank will absorb this energy and prevent any 
damage. 

Ultra capacitors are being considered over conventional capacitors due to their 
high capacitance relative to their small compact size. Since the technology is new, a short 
discussion on ultra capacitors will follow. 

Conventional capacitors store energy in the form of separated electric charge. The 
greater the area for storing charge and the closer the separated charges results in greater 
capacitance. A conventional capacitor achieves its cross-sectional area from plates of flat, 
conductive material. To realize high capacitance, this material can be wound in great 
lengths to increase its surface area. A conventional capacitor separates its charged plates 
with a dielectric material, sometimes plastic or paper film, or even a ceramic. These 
dielectrics can be made only as thin as the available films or other applied material. 

An ultra capacitor gets its surface area from a porous carbon-based electrode 
material. The porous structure of this material allows its surface area to approach 2000 
square meters per gram, much greater than can be accomplished using flat or textured 
films and plates. An ultra capacitor’s charge separation distance is determined by the size 
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of the ions in the electrolyte, which are attracted to the charged electrode. This charge 
separation is much smaller than can be accomplished using conventional dielectric 
materials. 

The combination of enormous surface area and extremely small charge separation 
gives the ultra capacitor its outstanding capacitance (1 Farad, 100 Volts) relative to 
conventional capacitors. 

E. Battery Background 

As shown in the candidate electrical system description, two battery banks, one 
connected to the Vetronics bus and the other to the auxiliary bus (see Figure 4-1), will be 
modeled and validated. Therefore, a brief discussion on basic battery concepts is 
included. 

A battery is a device that converts the chemical energy contained in its active 
materials directly into electric energy by means of an electrochemical oxidation-reduction 
(redox) reaction. This type of reaction involves the transfer of electrons from one material 
to another through an electric circuit. In a non-electrochemical redox reaction, such as 
rusting or burning, the transfer of electrons occurs directly and only heat is involved. 

While the term battery is often used, the basic electrochemical unit being referred 
to is the cell. A battery consists of one or more cells, connected in series or parallel, or 
both, depending on the desired output voltage or capacity. 

The cell consists of three major components. 
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1 . The anode or negative electrode - which gives up electrons to the external 
circuit and is oxidized during the electrochemical reaction. 

2. The cathode or positive electrode - which accepts electrons from the external 
circuit and is reduced during the electrochemical reaction. 

3. The electrolyte - the ionic conductor - which provides the medium for 
transfer of electrons, as ions, inside the cell between the anode and the 
cathode. The electrolyte is typically a liquid, such as water or other solvent, 
with dissolved salts, acids, or alkalis to impart ionic conductivity. 

The lead-acid battery, as will be used in the AAAV, uses lead dioxide as the 
active material for the positive electrode and metallic lead, in a high surface area porous 
structure, as the negative active material. The electrolyte is a sulfuric acid solution 
[Ref. 9], 
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V. 



COMPONENT MODELING AND VALIDATION 



A. Introduction 

This chapter contains a discussion of the modeling efforts applied to the NBC fan 
drive, ultra capacitor, and battery unit. It also provides results and conclusions for each 
model. 



B. NBC Fan Drive Model 

The circuit diagram of the NBC fan motor and its associated three-phase 
inverter drive is shown in Figure 5-1. The resistance and the inductance of the cable 
connecting the motor drive to the high power 28-volt bus are included in the diagram as 
nine and Z,iin e . The input capacitance of the inverter is given by Cdc- 




Figure 5-1, NBC Motor Drive System Circuit Diagram. 
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As discussed in Chapter 4, the equations describing the brushless dc motor are expressed 



as: 



v r qs = (r s + pL s ) i r qs + co r L/ ds + co r X' m 



(5.1) 



v r ds ={r s +pL s )i r ds -a) r Lj r qs 



(5.2) 



hV d\ 



T = 



\^y 









(5.3) 



where 




(5.4) 



and P is the number of poles. These equations together with the transformation of the 
phase voltages into the rotor reference frame and the inverse transformation of the rotor 
reference frame currents back to machine variables are coded in ACSL in a 
straightforward manner. 

The mechanical speed is established by integrating 






(5.5) 
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where T e is the electromagnetic torque, Ti is the load torque, and J i is the inertia of the 
motor and the connected fan load. It is assumed that the load torque is proportional to the 
square of the rotor speed. In particular, 

T l = K.al (5.6) 

The actual inverter is a 120° discontinuous-current inverter based on sensorless 
control. The current model being exercised by PCK&A is that of the NBC fan motor 
with a 1 80° continuous-current inverter using Hall-effect sensors to detect rotor position. 
The first step in developing an accurate model was to build a model similar to that being 
currently exercised. This model would allow for a better understanding of the inverter 
and motor assembly. 

The 180° inverter model was developed and exercised. The transistor gating 
signals are derived based on the waveforms illustrated in Figure 4-4. The phase voltages 
are then determined from (4.22) - (4.24). The current drawn by the inverter is 
established by summing the phase currents corresponding to the gated transistors 
connected to the top rail. The parameters used are shown in Table 5.1. The parameters 
were derived using manufacturer provided data. The complete ACSL model can be 
found in Appendix A. For convenience, a fixed 28-volt ideal source was used as the 
input to the inverter. 
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Symbol 


Value 


Description 


Origin 


4 


0.005 mH 


Stator self 
Inductance 


Calculated from manufacturer 
supplied data 

(average of upper and lower values) 


r s 


0.003 Q 


Stator resistance 


Calculated from manufacturer 
Supplied data 

(average of upper and lower values) 


r 

m 


0.00344 V-s/rad 


Back emf 
constant 


Calculated from manufacturer - 
supplied data 


P 


4 


Number of poles 


Communication from manufacturer 


C dc 


0.036 F 


Filter capacitance 


Communication from manufacturer 


^ line 


1 fjH 


Cable inductance 


Estimated 


r i,ne 


0.01 Q 


Cable resistance 


Estimated 


J m 


1.171X10“* kgm 2 


Inertia 


Estimated 


K m 


1.128X10' 6 Nms 


Load constant 


Estimated from manufacturer 
supplied data 



Table 5.1, NBC Fan Motor Drive Parameters. 



Figure 5-2 shows the first few cycles of the stator currents, /‘as, *bs, and ; cs , upon 
startup of the 180° inverter model using the data derived from communication with the 
manufacturer. With the low value of stator resistance, these currents illustrate a large 
startup transient, on the order of 4000A for /as- Since rated current for the drive is 
approximately 125 A, this current draw represents a significant transient load on the high- 
power bus. 
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Figure 5-2, Initial Stator Currents of NBC Fan Model 
with 180° Continuous-Current Inverter Upon Startup. 



Figure 5-3 shows simulation results for the drive extending from startup to steady 
state. The variables plotted include co rm - rotor speed in rads/s, T e - electromagnetic 

torque in Newton-meters, /'as - stator current in Amps, and i\\ - inverter input current in 
Amps. The large spike in inverter input current is readily apparent as are the considerable 
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current harmonics resulting from the applied six-step voltage waveforms. The motor 
speed settles out at a final value of about 1250 rad/sec. 
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Figure 5-3, Plot of Start-up of NBC Fan Model with 
180° Continuous-Current Inverter. 
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Figure 5-4 illustrates steady-state waveforms for the 1 80° inverter drive. In 
particular, the figure shows variables Vas - stator voltage in Volts, /as - stator current in 
Amps, and in - inverter current in Amps, T e - electromagnetic torque in Newton-meters, 

and 0) m — rotor speed in rads/s. The 6co r dominant harmonic in in is apparent and the 

peak value of /as appears to be on the order of 800A. This 800A amplitude is much 
larger than the anticipated 125 A, though the final speed of the device is rather close to the 
anticipated rated speed (1310 rad/sec). The difference may be due to the drive actually 
being a 120° inverter drive or to an inaccurate assumption on the motor parameters. 
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Figure 5-4, Plot of Steady State of NBC Fan Model with 
1 80° Continuous-Current Inverter. 
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Research then began on implementation of a 120° discontinuous-current inverter 
using a sensorless control scheme. As discussed in Chapter IV, paragraph C, subsection 
3.b, sensorless control can be initiated by indirectly assessing the rotor position using the 
back emf of the motor. The back emf can be monitored to detect the instant at which the 
back emf in the unexcited phase crosses zero or reaches a predetermined level. The 
instant at which this happens is then used to establish a switching sequence for 
commutation of the transistors in the inverter. 

The main nuance between the simulation of the 120° and the 180° inverter is 
establishing which devices are conducting. In particular, when neither transistor in the 
inverter leg is gated, then one of the diodes may be conducting or the phase could be 
open-circuited. In addition depending on the back emf, the open-circuited phase voltage 
may change and cause one of the diodes to turn on. The logic used in the simulation for 
establishing the status of Di and D 4 is illustrated below. 

IF ((Tl) .OR. (T4)) THEN 

Dl = .false. 

D4 = .false. 

ELSEIF (ias .GT. hold) THEN 

D4 = .true. 

Dl = .false. 

ELSEIF (ibs .LT. -hold) THEN 

Dl = .true. 

D4 = .false. 

ELSEIF ((Vagoc .LT. (Vdc + eps)) .OR. (Vagoc .GT. -eps)) THEN 

Dl = .false. 

D4 = .false. 
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ELSEIF (Vagoc .GT. (Vdc + eps)) THEN 



D1 = .true. 

ELSE 

D4 = .true. 

ENDIF 

Note, if neither transistor is gated and the phase current is positive, then D 4 must conduct. 
Conversely, if the current is negative, the Di must conduct. Otherwise, the back emf 
must be assessed. If it is greater than the input voltage (Vdc), then the upper diode is 
turned on. If it is less than zero volts, then D 4 is turned on. If the open-circuit voltage 
remains between zero volts and Vd C , then the phase remains open circuited. 

The complete ACSL code for the NBC Fan with the 120° discontinuous-current 
inverter can be found in Appendix B. The 120° inverter model was developed and 
exercised. The X' m parameter was adjusted through simulation to yield steady-state stator 
current waveforms consistent with manufacturer data. The complete set of parameters is 
documented in Table 5.2. Once again a fixed 28-volt source was used as the input to the 
inverter. 
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Symbol 


Value 


Description 


Origin 


4 


0.005 mH 


Stator self 
Inductance 


Calculated from manufacturer 
supplied data 

(average of upper and lower values) 


r s 


0.003 Q 


Stator resistance 


Calculated from manufacturer 
Supplied data 

(average of upper and lower values) 


K 

tn 


0.006 V-s/rad 


Back emf 
constant 


Adjusted by simulation 


P 


4 


Number of poles 


Communication from manufacturer 


c dc 


0.036 F 


Filter capacitance 


Communication from manufacturer 


L line 


1 fiH 


Cable inductance 


Estimated 


hne 


0.01 Q 


Cable resistance 


Estimated 


4 


1.171X10" 4 kg m 2 


Inertia 


Estimated 


K m 


1.1 28X10'^ ms 


Load constant 


Estimated from manufacturer 
supplied data 



Table 5.2, NBC Fan Motor Drive Parameters. 



Figure 5-5 shows the first few cycles of the stator currents, ias, /'b S , and i cs , upon 
startup of the 120° inverter model assuming that rotor position is available. Similar to 
the 1 80° model, the currents illustrate a large starting transient, reaching levels of 
approximately 2600A. Also, following the initial swing, the currents assume the classical 
120° waveform, where the open-circuit waveform is apparent. 
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Figure 5-5, NBC Fan Model Plot with 120° 
Discontinuous-Current Inverter. 



Figure 5-6 illustrates a complete startup to steady state and documents 
co m - rotor speed in rads/s, T e - electromagnetic torque in Newton-meters, - stator 
current in Amps, and in - inverter current in Amps. 
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Figure 5-6, Plot of Start-up of NBC Fan Model with 
120° Discontinuous-Current Inverter. 
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Figure 5-7 shows two cycles of the steady-state variables: - stator voltage in 

Volts, /as - stator current in Amps, and i\\ - inverter current in Amps, T e - 
electromagnetic torque in Newton-meters, and co m — rotor speed in rads/s. The inverter 

current once again exhibits a 6 co r fundamental variation, but it clearly has a different 

harmonic content than in the 1 80° model. The a-phase terminal voltage closely 
resembles the back emf as it should with the quasi -rectangular phase current waveforms. 
The a-phase stator current has a peak value of about 120A, inline with what the 
manufacturer data predicted. 
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Figure 5-7, Plot of the Steady State of NBC Fan Model with 
120° Discontinuous-Current Inverter. 



The simulation of the NBC fan model with 180° continuous-current inverter was 
run again using the parameters in Table 5.2 to illustrate the effect of varying X ' m . The 



62 





results can be found in Appendix C. Also, the simulation of the NBC fan model with 
120° discontinuous-current inverter was run again using the parameters in Table 5.1 to 
illustrate the effect of changing A' m . The results can be found in Appendix D. 

After successful implementation of the NBC fan model with 120° discontinuous- 
current inverter with a known rotor position, work began on implementing sensorless 
control of the motor. In order to accomplish this the open-circuit voltage signal would 
need to be processed into a signal that could be utilized to gate transistors Ti through T$ 
as discussed in Chapter IV, paragraph C, subsection 3d. Figure 5-8 depicts the signal 
processing block diagram that was implemented in order to process the open-circuit 
voltage. 




Figure 5-8, Open-Circuit Voltage Signal Processing Block Diagram. 
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where co c = 10,000 , a - 2500 , r = 0.004 . 

Gi, a low-pass filter, smoothes out the choppy terminal voltage waveform (vdiffc)- 
G 2 , an integrator, introduces a 90° phase shift. G 3 serves to block any dc component that 
may be introduced when the integrator starts. The zero-crossings of v co (together with 
analogous v ao and v bo signals from phase A and B) may be used for gating Tj-T, 5. Figure 
5-9 documents the start-up to steady-state waveforms from the block diagram with the 
machine operating with rotor position feedback. In Figure 5-9 the integrator is disabled 
until near the first peak of v^ c . This ensures that xp c does not have a large dc offset that 
would then need to settle out. Figure 5-10 documents the steady-state waveforms from 
the block diagram with the machine operating with rotor position feedback. The 90° 
phase shift between Vd,jf c (and Xfj c ) and xp c (and v c0 ) is readily apparent. 
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Figure 5-9, Plot of the Start-Up of NBC Fan Model with 
120° Discontinuous-Current Inverter. 
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The logical statements required to transition to sensorless operation are illustrated 

below. 



IF ((vco ,LT. 0.0) .AND. (vbo .LT. 0.0)) THEN 
Tla = . FALSE. 

T2a = .TRUE. 

T3a = .TRUE. 

T4a = .FALSE. 

T5a = .FALSE. 

T6a = .FALSE. 

ELSEIF ((vbo .GT. 0.0) .AND. (vao .GT. 0.0)) THEN 
Tla = .FALSE. 

T2a = .FALSE. 

T3a = .TRUE. 

T4a = .TRUE. 

T5a = .FALSE. 

T6a = .FALSE. 

ELSEIF ((vao .LT. 0.0) .AND. (vco .LT. 0.0)) THEN 
Tla = .FALSE. 

T2a = .FALSE. 

T3a = . FALSE. 

T4a = .TRUE. 

T5a = .TRUE. 

T6a = .FALSE. 

ELSEIF ((vco .GT. 0.0) .AND. (vbo .GT. 0.0)) THEN 
Tla = .FALSE. 

T2a = .FALSE. 

T3a = . FALSE. 

T4a = .FALSE. 

T5a = .TRUE. 

T6a = .TRUE. 

ELSEIF ((vbo .LT. 0.0) .AND. (vao .LT. 0.0)) THEN 
Tla = .TRUE. 

T2a = .FALSE. 

T3a = .FALSE. 

T4a = .FALSE. 

T5a = .FALSE. 

T6a = .TRUE. 

ELSEIF ((vao .GT. 0.0) .AND. (vco .GT. 0.0)) THEN 
Tla = .TRUE. 

T2a = .TRUE. 

T3a = .FALSE. 

T4a = .FALSE. 

T5a = .FALSE. 

T6a = .FALSE. 

ENDIF 

Details for the entire simulation are documented in Appendix B. 
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Figure 5-10, Steady-State Waveforms From Block Diagram with the 
Machine Operating with Rotor Position Feedback. 
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The following simulation study was conducted to verify the feasibility of 
sensorless operation. Rotor position information was used to start the machine (no back- 
emf voltages are initially available). At t=0.05, Figures 5-1 1 through 5-12 illustrate the 
signals derived from the actual rotor position information (Ti a through T6 a ). Clearly the 
signals Ti a through T^are very nearly in phase with Ti through Tg, respectively, and can 
be used to gate the transistors. Figure 5-14 depicts transitioning from using the rotor 
feedback signals to using the sensorless derived signals. The rotor maintains 
synchronism, with only a small drop in speed. 
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Figure 5-11, Plot of the Switching Signals of NBC Fan Model with 
120° Discontinuous-Current Inverter. 
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Figure 5-12, Plot of the Switching Signals of NBC Fan Model with 
120° Discontinuous-Current Inverter. 



70 



000Z fr0 : frfr : 80 ZZ 




o 
l o 



LO 

r>- 

LOo“ 



o 

o 



o 

LO 



LO 

r>- 

<E •. 
LO° 



o 

o 



o 

LO 



(Do' 



o 

o 



o 

LO 



CD L 



O 

o 



5.00000 



5.08333 



5. 16667 



5.25000 
* 10’ 2 



5.33333 



5.41667 



5.50000 



Figure 5-13, Plot of the Switching Signals of NBC Fan Model with 
120° Discontinuous-Current Inverter. 
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Figure 5-14, Plot of the Rotor Feedback Signals and 
the Sensorless Derived Signals. 

While the model introduced does implement sensorless control with good results, 
the problem of transitioning to sensorless control from rotor position feedback must still 
be resolved. This transition is an area that can be addressed in future wok. 
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C. Ultra Capacitor Model 

Pinnacle Research Institute, Inc. provided tabular data on their 100V, IF, 5000J 
Ultra-Capacitor. The manufacturer’s measured data was entered into a MATLAB file, 
ultracap.m. This file is included as Appendix C. The following plot illustrates the 
manufacturer’s measured impedance versus frequency data. 



Impedance vs. Frequency 




Figure 5-15, Ultra Capacitor Impedance vs. Frequency 
(Manufacturer’s Measured Data). 

A computer model of the ultra capacitor would need to produce comparable 
results when exercised over the same frequency range. An attempt to duplicate the actual 
results would be initiated using an equivalent series resistance (ESR) as in Figure 5-16. 



73 



R 



eq 






1/SC 



Figure 5-16, ESR Circuit Diagram. 



The input impedance of the circuit is found by inspection to be 



which can be rewritten as 



eq 




o 

SCj 









w 



R\S + 



RC 



(5.10) 



(5.11) 



J 



An initial guess for the constants R and 1/RC was formulated using the PRI, Inc. 
frequency response data. After many iterations of manipulating these constants the 
following plot of impedance vs. frequency was achieved. 
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Figure 5-17, Ultra Capacitor Impedance vs. Frequency (Model). 



Figure 5-18 illustrates both actual data points and model generated data and 
confirms that the model is valid as a reasonable representation of the actual ultra 
capacitor. 
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Figure 5-18, Ultra Capacitor Impedance vs. Frequency (Actual and Model). 

The constants R and — were found to be: 

RC 

R= 0.045 ohms (5.12) 

— =5.35 s' 1 (5.13) 

RC 



and the transfer function representation for the ultra capacitor would be characterized as 



Z 



eq 



0.045(5 + 5.35) 

S 



(5.14) 
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D. Battery Model 



The AAAV electrical system consists of two battery banks, one connected to the 
Vetronics bus and the other to the Auxiliary bus. The bank connected to the Vetronics 
bus consists of two series connected 12- Volt (nominal) batteries. The bank connected to 
the auxiliary bus consists of three parallel branches with two series-connected 12 - Volt 
batteries per branch. In the current Purdue Distributed Computing AAAV system 
representation, the batteries are modeled as a Thevenin equivalent circuit, where the 
open-circuit voltage and the Thevenin resistance can be adjusted to correspond to the 
chemical and thermodynamic state of the battery. The resulting circuit diagram is shown 
in Figure 5-19. The associated parameters are summarized in Table 5.3. In Figure 5-19, 
Si, S 2 , and S 3 are closed for starting. After start-up S 2 and S 3 are opened, and Si is left 
closed during operation. 

When no battery is connected to the Vetronics bus, a bus voltage oscillation was 
observed. The stiffening action of the battery seems to mitigate this instability. It is 
hypothesized that the oscillation is due to a limit cycle imposed by the generator voltage 
regulator, particularly at low power. 
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Figure 5-19, Circuit Diagram of Battery Banks. 



Symbol 


Description 


Value 


Origin 


R&; 


Vetronics bus battery series 
resistance 


0.007648 Q 


Calculated using 
charge/discharge data 


P-62 


Auxiliary bus battery series 
resistance 


0.00255 Q 


Calculated using 
charge/discharge data 


R 12 


Resistance of cable between 
main and auxiliary buses 


0.0002 Q 


Estimated 


Lbi 


Vetronics bus battery series 
inductance 


1.0 pH 


Estimated 


1-62 


Auxiliary bus battery series 
inductance 


1.0 pH 


Estimated 



Table 5.3, Battery Parameters. 
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The model of Figure 5.19 is preliminary and further laboratory testing has been 
conducted on the 800U Optima battery. Of particular interest was the battery’s frequency 
response at various levels of charge. The circuit shown in Figure 5.20 was constructed 
and the parameters for the test circuit are summarized in Table 5.4. The purpose of Ri is 
to establish a quiescent load condition on the battery. Transistor Qi is gated at the user- 
specified frequency, with a 50% duty cycle. With R .2 larger than Ri, switching Qi results 
in a smaller rectangular sink of current from the battery. Therefore, i te rm consists of a 
steady-state dc level with a small rectangular ripple. This AC component of current flows 
through the battery resulting in an AC component of voltage across the battery (v te rm)- 
Therefore, if the assumed internal battery model is that given by Figure 5-21, the AC 
component of the measured voltage will facilitate the calculation of the impedance given 
by rbatt and Lbatt- The current was monitored with a Tektronix, 8-Amp, clamp-on, current 
probe, while the terminal voltage was monitored with a direct connection to the 
oscilloscope and with a digital handheld multi-meter. Both signals were displayed on a 
Tektronix, 60 MHz, digital storage oscilloscope. The ratio of the amplitudes of the 
ripples in v tenn and i te rm provides a measure of the impedance magnitude, while the phase 
angle (0) between those signals relates how the impedance divides between real and 
reactive components. Figure 5-21 reflects the assumed internal battery representation. 
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Figure 5.20, Circuit Diagram for Optima Battery Frequency Response Test. 
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Figure 5-21, Assumed Internal Battery Representation. 
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Symbol 


Description 


Measured 

Value 


General Information 


R. 


Ceramic Resistor 


4.0 Q 


4 X 1 6 Q resistors in parallel, 
1 0 W each 


r 2 


Carbon Resistor 


40.1 Q 


(1) 39 Q 


Qi 


Power MOSFET (2N6796) 




8.0 A, 100 V 




Optima 800U Battery 


12.8 


56 A-hr 



Table 5.4, Test Circuit Parameters. 

Testing was conducted at frequencies of 20 Hz to 20 kHz. The upper end of 
20kHz was deemed sufficient since the ripple frequency out of the 1 2-pulse 
generator/rectifier is less than 20 kHz. The testing was conducted at three different charge 
states, (100%, 75%, and 50%). The charge state of the battery is determined using the 
manufacturer supplied capacity of 56 A-hr. Full battery specifications are provided in the 
Manufacturer Data Sheet listed in Appendix F. Given that 

\A-hr-=\—X 3600 5 = 3600 C , (5.15) 

s 

the stored charge (q s tored) is determined as 



q slored =56 A-hr X 3600 — — = 201,600 C . 

A-hr 



(5.16) 
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Since charge is related to current by 



?={'<*■ (5-17) 

then for constant current ( 1 ), equation (5-17) becomes 

q = I At (5.18) 

By drawing a fixed current (I) over a specific time interval (At), the amount of charge left 
in the battery can be estimated. 

Initial test results of the battery at full charge, V = 12.8 Volts, I =3 .5 A, and / = 
2kHz are depicted in Figure 5-22. 
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Figure 5-22, Test Results with Battery at Full Charge, 

V = 12.8 Volts, I =3. 5 A, and/= 2kHz 
CHI = current (1 Amp/div), CH2 = voltage (lVolt/div). 

The plot reveals V (mV pp ) - 3.0, I (mA pp ) = 3 10, 0 (rads), = 7 t, giving an internal 
resistance, R (mQ) = 9.67, and internal inductance, L (mH) = 0.0. 

The battery was then drained at a current level of approximately 15.5 A for 52 
minutes to arrive at a charge state of about 75% rated. At 75% charge test results of the 
battery with V = 1 1 .9 Volts, I =3.3A, and/= 2kHz are depicted in Figure 5-23. 
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Figure 5-23, Test Results with Battery at 75% Charge, 
V = 1 1.9 Volts, I =3. 3 A, and / = 2kHz 
CHI = current (1 Amp/div), CH2 = voltage (lVolt/div). 



The plot reveals V (mV pp ) = 2.5, / (mA pp ) = 300, 0 (rads), = it, giving an internal 
resistance, R (mQ) = 8.33, and internal inductance, L (mH) = 0.0. 

The battery was then drained at a current level of approximately 15A for 55 
minutes to arrive at a charge state of about 50% rated. At 50% charge test results of the 
battery with, V = 1 1.8 Volts, I =3.25 A, and /= 2kHz are depicted in Figure 5-24. 
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Figure 5-24, Test Results with Battery at 50% Charge, 
V = 1 1.8 Volts, I =3 .25 A, and/= 2kHz 
CHI = current (1 Amp/div), CH2 = voltage (1 Volt/div). 



The plot reveals V (mV pp ) = 2.7, / (mA pp ) = 290, 0 (rads), = n, giving an internal 
resistance, R (mQ) = 9.31 , and internal inductance, L (mH) = 0.0. 

Complete laboratory test results for all frequencies are summarized in Table 5.5. 
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100% Charge 


75% Charge 


50% Charge 


/ 


V 

mVpp 


i 

mApp 


0 

rad 


R 

mft 


L 

mH 


V 

mVpp 


I 

mApp 


e 

rad 


R 

mH 


L 

mH 


V 

mVpp 


i 

mApp 


0 

rad 


R 

mfl 


L 

mH 


20 

Hz 


3.1 


310 


n 


10.0 


0.0 


2.5 


260 


7t 


9.62 


0.0 


2.8 


280 


71 


10.0 


0.0 


200 

Hz 


2.9 


300 


n 


9.67 


0.0 


2.5 


300 


71 


8.33 


0.0 


2.7 


280 


71 


9.64 


0.0 


2 

kHz 


3.0 


310 


n 


9.67 


0.0 


2.5 


300 


71 


8.33 


0.0 


2.7 


290 


71 


9.31 


0.0 


20 

kHz 


3.5 


370 


n 


9.46 


0.0 


2.8 


320 


7t 


8.75 


0.0 


2.9 


310 


71 


9.35 


0.0 



Table 5.5, Laboratory Test Results. 



Test results conclude that the internal resistance of the battery does exhibit a small 
change as a result of frequency and charge-state of the battery. Further, tests conclude 
that the internal inductance of the battery is negligibly small. It appears reasonable that 
the battery can be adequately represented by an equation of the form 

V /erm = V i„. “V/e™ (5-19) 



If v lerm is the same as the bus voltage and it is a state variable then 



l term 



v — V, 

int bus 



'bait 



(5-20) 



where i lerm is then available to evaluate the derivative of the bus voltage. 



86 



VI. CONCLUSION 



A. Synopsis of Work Completed 

1. Introduction 

The Advanced Amphibious Assault Vehicle (AAAV) was introduced 
and a general system description was set forth. The introduction and description are 
critical to understanding where and how the AAAV fits into the Marine Corps 
amphibious triad (AAAV, V-22, LCAC) to execute the concepts of Operational 
Maneuver from the Sea and Ship to Objective Maneuver. The AAAV’s operating modes 
and associated mode characteristics were also documented and developed. 

2. Electrical System Architecture 

The entire electrical system architecture is introduced in block diagram 
form and the system’s primary components are listed and discussed. A discussion on the 
theory behind distributed computing techniques through networking is included. This 
technique is the foundation of P. C. Krause and Associates’ troubleshooting tool that can 
be used by the Marine Corps in addressing current electrical system issues that must be 
resolved. An Advanced Continuous Simulation Language (ACSL) background is ' 
included. This software is used in this research to model and simulate various 
components of the electrical system. 
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3. 



Candidate Electrical System Architecture 



From the entire electrical system architecture, a scaled down candidate 
system architecture is developed. Components from the candidate system are tested, 
modeled, simulated, and validated. The Nuclear, Biological, and Chemical (NBC) 
fan/filter motor drive, an ultra capacitor, and a battery unit are the components selected 
for research in this document. 

4. NBC Fan/Filter Motor Drive 

A simulation of the NBC fan motor drive with a 1 80° continuous- 
current inverter using Hall-effect sensors to detect rotor position was modeled in ACSL. 
The simulation allowed for a thorough understanding of the inverter and motor assembly, 
and also served as a tool for gaining an understanding of ACSL simulations and 
representative signal waveforms. 

Efforts then shifted to implementing a NBC fan motor drive model with the 120° 
discontinuous-current inverter based on sensorless control. This model was achieved in 
two stages. First, the 120° discontinuous-current inverter model was developed 
assuming the availability of rotor position. The second step was to incorporate sensorless 
control into the simulation. This was achieved by indirectly sensing the rotor position 
using the back emf of the motor. The switching sequence for commutation of the 
transistors in the inverter was established using this back emf. The results were 
consistent with theoretical expectations ensuring this model was a valid representation of 
the NBC fan motor and inverter assembly. 
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5. 



Ultra Capacitor 



Tabular data was collected from Pinnacle Research Institute Inc., on 
the 100V, IF, 5000J, ultra capacitor. This data was analyzed and plotted using MATLAB 
and an impedance vs. frequency plot was developed. An attempt to duplicate this plot 
using a model with an equivalent series resistance (ESR) was initiated. A transfer 
function was developed and an initial guess was made for constants. By manipulating the 
transfer function constants, an impedance vs. frequency plot matching that plot from the 
manufacturer’s data was achieved. 

6. Battery Unit 

The area of concern for the battery unit was its frequency response 
characteristics at different charge states. A test circuit was constructed that would load 
the battery drawing a constant DC level of current and an AC component of current. The 
current and the voltage levels were displayed on an oscilloscope and monitored. From 
the steady-state and ripple current and voltage waveforms, an internal impedance was 
calculated. These tests were conducted at three charge states (100%, 75%, and 50%) of 
the battery over a frequency range of 20 Hz to 20 kHz. The upper frequency level of 20 
kHz was based on the ripple frequency of the 12-pulse generator/rectifier. Test results 
concluded that the internal impedance of the battery was almost purely resistive with a 
negligibly small internal inductance. Further, tests concluded that the internal resistance 
of the battery did fluctuate as a result of charge state and frequency. 
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B. DOD Relevance 



The AAAV will join the Marine Corps inventory to complete the amphibious 
triad which also includes the V-22 Osprey, tilt-rotor aircraft, and the Landing Craft, Air 
Cushion (LCAC). The AAAV will provide the principal means of water mobility, land 
mobility, and direct fire support to Marine Corps infantry units. 

As the AAAV travels through its prototyping stage in an effort to successfully 
pass its operational testing, electrical system issues and problems have arisen which must 
be addressed before entering into full-scale production. Computer simulation of the 
electrical system will facilitate the analysis and evaluation of hardware design prior to 
and in concert with actual production. Accurate simulation will provide invaluable 
insight into design flaws or system limitations before extensive capital investments are 
made on hardware. 

C. Further Work 

Complete simulation for the AAAV’s entire electrical system is the next step in 
providing the Marine Corps with a tool where by it can successfully analyze and evaluate 
the system. Once a complete and accurate representative model exists, a distributed 
computing network must be established in order to exercise the entire model. Funding 
for such a network at NPS has not yet been obtained and a future researcher may be able 
to pursue possible avenues for funding. 

Other research that can begin now includes a paralleling strategy for the auxiliary 
generator. The current auxiliary generator is incapable of supporting predicted power 
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requirements. Finally, research could be initiated in an attempt to resolve instability 
problems in the electrical system during low power operation. General Dynamic’s 
engineers recently discovered low power instability problems during prototype testing. 
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APPENDIX A. [ACSL CODE FOR NBC FAN MODEL WITH 180° INVERTER] 



PROGRAM NBC Fan Model w/ 180 Degree Inverter 
INITIAL 



MAXTERVAL rnaxt = 1.0e-6 
MINTERVAL mint = 1.0e-7 
CINTERVAL cint = 1.0e-4 
ALGORITHM ialg = 5 
CONSTANT "ts top = 0.165 



• initial circuit constants • 



CONSTANT rL=0 . 01 
CONSTANT LL=. 000001 
CONSTANT Cdc=. 036 
CONSTANT Vdc=28 . 0 

LOGICAL SA 
LOGICAL S3 
LOGICAL SC 



" switch initial conditions * 



SA = .false. 
S3 = .false. 
SC = .false. 



END ! "of initial" 



DYNAMIC 



TERMT Ct .GE. tstop) 



DERIVATIVE 



{"Switch Status Generation" 

CONSTANT ?i=3 . 14159265 
twopi = 2.0*pi 
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PROCEDURAL ( the tar = thetarm) 



thetar = (Poles/2.) * thetarm 

SCHEDULE reset .XP. ( thetar- twopi ) 
SCHEDULE reset2 .XN. (thetar) 



END ! End of Procedural 



["Firing of Switches" 



IF ((thetar .GE. 0.0) .AND. (thetar .LT. pi/6.)) THEN 
SA = . true . 

SB = .false. 

SC = . false . 

ELSEIF ((thetar .GE. pi/6.) .AND. (thetar .LT. pi/2.)) THEN 
SA = .true. 

SB = . true . 

SC = . false . 



ELSEIF ((thetar .GE. (pi/2.)) .AND. (thetar .LT. (5.*pi/6.))) 
SA = . false. 

SB = .true. 

SC = .false. 

ELSEIF ((thetar .GE. (5.*pi/6.)) .AND. (thetar .LT. (7.*pi/6. 
SA = .false. 

SB = .true. 

• SC = . true . 

ELSEIF ((thetar .GE. (7.*pi/6.)) .AND. (thetar .LT. (3.*pi/2. 
SA = .false. 

SB = .false. 

SC = .true. 



ELSEIF ((thetar .GE. (3.*pi/2.)) .AND. (thetar .LT. (ll.*pi/S. 
SA = . true . 

S3 = .false. 

SC - . true . 



ELSEIF ((thetar .GE. (ll.*?i/6.)) .AND. ( the car .LT. (2.*pi))) 
SA = . true . 

SB = .false. 

SC = .false. 

ENDIF 



! "Assume lower switches are conducting" 



THEN 



) ) THEN 



) ) THEN 



) ) ) THEN 



THEN 
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PROCEDURAL (Vas , Vbs, Vcs, iii = SA, SB, SC, Vdc, ias, ibs, ics) 



Vap = 0.0 
iia = 0.0 
Vbp =0.0 
iib = 0.0 
Vcp = 0.0 
iic = 0.0 

! " Check if Cop switch is closed" 

IF (SA) THEN 

Vap = Vdc 
iia = ias 

ENDIF 

IF (SB) THEN 

Vbp = Vdc 
iib = ibs 

ENDIF 

IF (SC) THEN 

Vcp = Vdc 
iic = ics 

ENDIF 

!" Establish the phase voltages" 

Vnp = (Vap + Vbp + Vcp) /3.0 

Vas = Vap - Vnp 
Vbs = Vbp - Vnp 
Vcs = Vcp - Vnp 

!" Find the inverter input current" 
iii = iia + iib + iic 
END ! End of Procedural 



I “ Link Dynamic Equations " 
Vcap =28.0 



! "Brushless DC Motor Constants'* 
Ls = .000005 
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r s = 0.003 
lamref = 0.006 
Poles = 4.0 
Jm = 0.00011711 
Km = 0.000001128 

! "Brushless DC Motor Equations* 

'"Rotor Ref Frame Voltages" 

Vqsr = 2 . *cos ( thetar) *Vas/3 . + 2 . *cos ( thetar-2 . *pi/3 . ) *Vbs/3 . & 
+ 2 . *cos (thetar + 2 . *pi/3 . ) *Vcs/3 . 



Vdsr = 2 . *sin ( thetar) *Vas/3 . +■ 2 . *sin( thetar-2 . *pi/3 . ) *Vbs/3 . & 
+ 2 . *sin( thetar + 2 . *pi/3 . ) *Vcs/3 . 



! " Inverter Currents " 

ias = cos (thetar) *iqsr + sin( thetar) *idsr 

ibs = cos (thetar - 2 . *pi/3 . ) *iqsr + sin(thetar - 2 . *pi/3 . ) *idsr 
ics= -ias-ibs 



!• Rotor Ref Frame Currents" 

piqsr = (-rs*iqsr - wr*Ls*idsr - wr*lamref + vqsr) /Ls 
iqsr = INTEG (piqsr , 0.0) 

pidsr = ( -rs*idsr + wr *Ls * iqsr + vdsr) / Ls 
idsr = INTEG (pidsr ,0.0) 

'"Rotor Ref Frame Angle Thetar and wr" 
pwrm - (Te - T1 ) / Jm 
wrm = INTEG (pwrm, 0.0) 
wr = (Poles/2.) * wrm 
p the t arm = wrm 

thetarm = INTEG (pthetarm , 0.0) 



! " Torque Equations" 
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Te = (Poles/2.) * ( (3 . /2 . ) # lainref) * iqsr 
T1 = Km * (wrm*wrm) 

END ! " End of derivative " 



DISCRETE reset 

the t arm = 0.0 

END 

DISCRETE *reset2 

the t arm = the t arm + twopi 

END 



END ! " End of dynamic section" 
END ! "End of program" 
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APPENDIX B. [ACSL CODE FOR NBC FAN MODEL WITH 120° INVERTER] 



PROGRAM NCB FAN MODEL (AAAV) 
INITIAL 



{"Preliminary stuff" 
MAXTERVAL maxt = 1.0e-7 
MINTERVAL mint = 1.0e-7 
CINTERVAL cint = 1.0e-5 
ALGORITHM ialg = 5 
CONSTANT tstop = 0.007 



! "max time step" 

! "min time step" 
‘"communication interval" 
! "4th order Runge-Kutta" 

! "adjustable stop- time" 



** initial circuit constants-- 

CONSTANT. rL=0 . 01 
CONSTANT LL=. 000001 
CONSTANT Cdc= .036 
CONSTANT Vdc=28 . 0 

! " parameters for determining 

CONSTANT hold = l.e-1 
CONSTANT eps = 0.7 

!"Brushless DC Motor Constants" 
CONSTANT Ls = .000005 
CONSTANT rs = 0.003 
CONSTANT lamref = 0.006. 

CONSTANT Poles =4.0 
CONSTANT Jm = 0.00011711 
CONSTANT Km = 0.000001128 



{"input line series R" 

{"input line series L" 

{"input capacitance" 

{"input voltage" 

circuit status" 

{"holding current for diodes" 
{ “diode drop" 



{"self inductance" 
{"stator resistance" 
i "magnet parameter" 
{"number of poles" 

{ "motor+load inertia" 
! "TL=Km*wrm*wrm" 



" logical variables for transistor conduction status" 

LOGICAL Tl , T2 , T3 , T4 , T5 , T6 

" logical variables for diode conduction status" 

LOGICAL Dl, D2, D3, D4, D5, D6 



" logical variables derived from sensorless signals" 

LOGICAL Tla, T2a, T3a, T4a, TSa, TSa 

" switch initial conditions " 

" all devices off" 

Tl = . false . 

T2 = . false. 

T3 = . false . 

T4 = . false . 

T5 = .false. 

T6 = .false. 

Dl = . false . 

D2 = .false. 

D3 = .false. 

D4 = . false . 

D5 = .false. 

D6 - . false . 
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all sensorless signals initially false" 

Tla = .false. 

T2a = .false. 

T3a = .false. 

T4a = .false. 

T5a = .false. 

T6a = .false. 

" logical signals initiate integration of phase voltages 

LOGICAL flaga, flagb, flagc 
flaga = .TRUE, 
flagb = .TRUE, 
flagc = .TRUE. 

" logical symbol to switch over to sensorless signals" 

LOGICAL NOSENSE 
NOSENSE = .TRUE. 

END ! “of initial" 

DYNAMIC 



• simulation termination condition" 

TERMT (t .GE. tstop) 

DERIVATIVE 

‘"constants needed for manipulating angles" 
CONSTANT pi=3 .14159265 
twopi = 2.0*pi 

CONSTANT alpha =1.047198 !"60 degrees" 

PROCEDURAL ( the tar = thetarm) 

!" — shift rotor electrical angle by alpha" 
thetarx= (Poles/2.) * thetarm + alpha 

! "keep angle between 0 and twopi" 

IF (thetarx .GE. twopi) THEN 

thetarx= thetarx - twopi 
ELSEIF (thetarx . LT. 0.0) THEN 

thetarx = thetarx + twopi 

ENDIF 

the tarx= MOD (thetarx, twopi) 

! "recover the tar" 
thetar = thetarx - alpha 
END ! " End of Procedural" 



! "Firing of Switches" 

!"if NOSENSE is TRUE, use rotor position signals for Tl-TS 
! " if NO SENSE is FALSE, use open-circuit-derived signals" 

IF (NOSENSE) THEN 

IF ((thetarx .GE. 0.0) .AND. (thetarx .LT. pi/ 3.)) THEN 

T1 = .true. 

T2 = .false. 

T3 = .false. 
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(thetarx . LT. (2 . *pi/3 . ) ) ) THEN 



T4 = .false. 

T5 = .false. 

T6 = . true . 

ELSEIF ((thetarx .GE. pi/3.) .AND. 

T1 = .true. 

T2 = . true . 

T3 = .false. 

T4 = . fa’lse. 

T5 = .false. 

T6 = .false. 

ELSEIF ((thetarx .GE. (2.*pi/3.)) .AND. (thetarx .LT. pi 
T1 = .false. 

T2 = . true . 

T3 = .true. 

T4 = .false. 

T5 = .false. 

T6 = .false. 

ELSEIF ((thetarx .GE. pi) .AND. (thetarx .LT. (4.*pi/3.)) 
T1 = .false. 

T2 = .false. 

T3 = . true . 

T4 = .true. 

T5 = .false. 

T6 = .false. 



ELSEIF ((thetarx .GE. (4. ♦pi/3.)) .AND. (thetarx .LT. (5. 
T1 = .false. 

T2 = .false. 

T3 = .false. 

T4 •= . true . 

T5 = .true. 

T6 = .false. 

ELSEIF ((thetarx .GE. (5. ♦pi/3.).) .AND. (thetarx .LT. (2. 
T1 = .false. 

T2 = .false. 

T3 = .false. 

T4 = .false. 

T5 = .'true. 

T6 = . true . 

ENDIF 

ELSE 

T1 = Tla 
T2 = T2a 
T3 = T3a 
T4 = T4a 
T5 = T5a 
T6 = T5a 

ENDIF 



! " Determine diode conduction status" 

I " establish the open-circuit voltage " 

Vagoc = ( ( Vdc / 2 . ) - ( ( 3 . ♦lamref ♦wr’cos ( thetar ) ) 12 . 
Vbgoc = ((Vdc/2.) ^ (( 3 . ♦lamref *wr*cos ( thetar - (2 
Vcgoc = ((Vdc/2.) + ({ 3 . ♦lamref w wr*cos ( thetar + (2 



) ) THEN 



) THEN 



♦pi/3.))) THEN 



♦pi))) THEN 



) ) 

. ♦pi/3 . ) ) ) / 2 . ) ) 
. ♦? i / 3 . ) ) )/2.) ) 
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i "diode logic for PHASE LEG A” 

l"if either transistor is gated, then either the transistor " 

“ or diode will be conducting" 

IF ((Tl) .OR. (T4 ) ) THEN 
D1 = . false . . 

D4 = . false . 

!"if the leg transistors are not gated:" 

!"if ias>hold, then the lower diode is on, upper is off" 

ELSEIF (ias .GT. hold) THEN 
D4 = . true . 

D1 = .false. 

i"if ias<-hold, then the upper diode is on, lower is off" 

ELSEIF (ias . LT. -hold) THEN 
D1 = . true . 

D4 = .false. 

!"if open circuit voltages are between Vdc + eps and -eps" 

!"then both diodes are off" 

ELSEIF ( (Vagoc .LT. (Vdc + eps)) .OR. (Vagoc .GT. -eps)) THEN 
D1 = .false. 

D4 = .false. 

!"if open circuit voltage is above vdc + eps, upper diode is made" 
!"to turn on" 

ELSEIF (Vagoc .GT. (Vdc + eps)) THEN 
D1 = . true . 

! "otherwise the lower diode turns on" 

ELSE 

D4 = . true . 

ENDIF 



'"diode logic for PHASE LEG B" 

IF ( (T3 ) -OR. (T6) ) THEN 
D3 = '.false. 

D6 = .false. 

ELSEIF (ibs .GT. hold) THEN 
D3 = .false. 

D6 = .true. 

ELSEIF (ibs . LT . -hold) THEN 
D3 = . true . 

Do = .false. 

ELSEIF ( (Vbgoc .LT. (Vdc + eps)) .OR. (Vbgoc .GT. -eps)) THEN 
D3 = .false. 

D6 = .false. 

ELSEIF (Vbgoc .GT. (Vdc + eps)) THEN 
D3 = .true. 



102 



ELSE 

D6 = .true. 

ENDIF 



! "diode logic for PHASE LEG C" 

IF' ( (T5) .OR. (T2 ) ) THEN 
D5 = .false. 

D2 = .false. 

ELSEIF ( ics .GT. hold) THEN 
D5 = .false. 

D2 = . true . 

ELSEIF (ics .LT. -hold) THEN 
D2 = .false. 

D5 = .true. 

ELSEIF ( (Vcgoc .LT. (Vdc + eps ) ) .OR. (Vcgoc .GT. -eps)) THEN 
D5 = .false. 

D2 = -false. 

ELSEIF (Vcgoc .GT. (Vdc + eps)) THEN 
D5 = .true. 

ELSE 

D2 = . true . 

ENDIF 



’"Procedural for finding stator voltages" 

PROCEDURAL (Vas, Vbs , Vcs , iii = Tl, T2 , T3 , T4, T5 , T6 , Dl, D2 , D3 , D4 , 
D5, D6, ias, ibs , ics) 



! " Finding pole voltages and currents for A-phase" 

! "upper device is on" 

IF \ (Tl) .OR. (Dl) ) THEN 
Vag = Vdc 
iia = ias 

!" lower device is on" 

ELSEIF ( (T4 ) .OR. (D4 ) ) THEN 
Vag =0.0 
iia = 0.0 

i " leg is open circuited" 

ELSE 

Vag = ((Vdc/2.) + ( (3 . * lamref *wr *cos ( the tar) ) /2 . ) ) 
iia = 0.0 

ENDIF 
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! " Finding pole voltages and currents for B-phase" 

IF ( (T3 ) . OR. (D3 ) ) THEN 

Vbg = Vdc 
iib = ibs 

ELSEIF ( (T6 ) .OR. ( D6 ) ) THEN 
Vbg =0.0 
iib = 0.0 

ELSE 

Vbg = ((Vdc/2.) + { { 3 . *lamref *wr *cos ( thetar - (2 . *pi/3 . ) ) ) 12 . ) ) 
iib = 0.0 

ENDIF 



! " Finding pole voltages and currents for C-phase" 

IF ( (T5) .OR. (D5 ) ) THEN 
Vcg = Vdc 
iic = ics 

ELSEIF ( (T2 ) . OR. ( D2 ) ) THEN 

Vcg =0.0 
iic = 0.0 

ELSE 

Vcg = ((Vdc/2.) + ( (3 . *lamref *wr*cos (thetar + (2 . *pi/3 . ) ) ) /2 . ) ) 
iic = 0.0 

ENDIF 



! " Establish the phase voltages" 

Vas = { ( 2 . *Vag/3 . ) - (Vbg/3.) - (Vcg/3.)) 

Vbs = ( ( 2 . *Vbg/3 . ) - (Vag/3.) - (Vcg/3.)) 

Vcs = ( ( 2 . *Vcg/ 3 . ) - (Vag/3.) - (Vbg/3.)) 

1 " Find the inverter input current" 

iii = iia + iib + iic 

END ! " End of Procedural" 

! "evaluate the neutral voltage" 
vng = (vag + vbg + vcg)/ 3.0 



! "C-Phase Sensorless Signal derivation" 

CONSTANT wc = 10000.0 !" low-pass filter frequency" 

CONSTANT tau2 = 0.004 

f "calculate the terminal voltage: should be Vcs" 
vdiffc = vcg - vng 

!"lowpass filter to eliminate abrupt transitions" 

CONSTANT xflcic =0.0 
pxflc = -wc*xflc + wc*vdiffc 
xflc = INTEG (pxflc , xflcic) 
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! "activate integrator following first peak" 

SCHEDULE intone .XP. (16.0 - xflc) 

IF (flagc) then 

pxf2c = 0.0 

ELSE 

pxf2c = 2500.0*xflc 

ENDIF 

xf 2c = INTEG (pxf 2c , 0.0) 

! "apply result to third filter to eliminate dc component 
pvex = (xf2c - vcx)/tau2 
vex = INTEG (pvex, 0.0) 
vco = -vex + xf2c 

! "do same for the A-Phase" 
vdiffa = vag - vng 
CONSTANT xf laic =0.0 
pxfla = -wc*xfla + wc*vdiffa 
xfla = INTEG (pxf la, xflaic) 

SCHEDULE intona .XP. (16.0 - xfla) 

IF (flaga) then 

pxf 2a = 0.0 

ELSE 

pxf 2a = 2500 . 0*xfla 

ENDIF 

xf 2a = INTEG (pxf 2a, 0.0) 
pvax = (xf2a - vax)/tau2 
vax = INTEG ( pvax , 0.0) 
vao = -vax + xf2a 

! "do- same for the B-Phase" 
vdiffb - vbg - vng 
CONSTANT xflbic =0.0 
pxflb = -wc*xflb + wc*vdiffb 
xflb = INTEG (pxflb, xflbic) 

SCHEDULE intonb .XP. (16.0 - xflb) 

IF (flagb) then 

pxf 2b = 0.0 

ELSE 

pxf2b = 2500 . 0*xf lb 

ENDIF 

xf 2b = INTEG (pxf 2b, 0.0) 
pvbx = (xf2b - vbx)/tau2 
vbx = INTEG (pvbx, 0.0) 
vbo = -vbx + xf2b 



vcsoc = lamref *wr*cos ( thetar + 2 . *pi/ 3 . ) 



! " Link Dynamic Equations" 

{"removed for simulation development convenience" 
IpiL = (Vac - rL*iL - Vcap) /LL 
! iL = INTEG (piL , 0.0) 

IpVcap = (iL - iii)/Cdc 
IVcap = INTEG (pVcap , 23.0) 

Vcap = 28.0 ! "assume hard input voltage 
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! "Brushless DC Motor Equations” 

’"Rotor Ref Frame Voltages” 

Vqsr = 2 . *cos { the tar ) *Vas/3 . + 2 . *cos ( thetar-2 . *pi/3 . ) *Vbs/3 . & 
+ 2 . *cos { thetar + 2 . *pi/3 . ) *Vcs/3 . 

Vdsr = 2 . *sin( thetar) *Vas/3 . + 2 . *sin ( thetar-2 . *pi/3 .) *Vbs/3 . & 
+ 2 . *sin ( thetar + 2 . *pi/3 . ) *Vcs/3 . 

! ” Phase Currents ” 

ias = cos (thetar) *iqsr + sin ( thetar) *idsr 

ibs = cos (thetar - 2 . *pi/3 . ) *iqsr + sin (thetar - 2 . *pi/3 . ) *idsr 
ics= -ias-ibs 



! "Differential equation for the machine” 

piqsr ={-rs*iqsr - wr*Ls*idsr - wr*lamref + vqsr)/Ls 

iqsr = INTEG (piqsr , 0.0) 

pidsr = {-rs*idsr + wr*Ls*iqsr + vdsr) /Ls 
idsr = INTEG (pidsr , 0.0) 

!" Torque and mechanical dynamics" 

Te = (Poles/2.) * ((3./2.) * lamref) * iqsr 
Tl = Km * (wrm* wrm) 

pwrm = (Te - Tl)/ Jm 
wrm = INTEG (pwrm, 0.0) 
wr = (Poles/2.) * wrm 



pthetarm = wrm 

thetarm = INTEG (pthetarm , 0.0) 

! "use sensorless signals to derive switching signals” 
IF ( { vco .LT. 0.0) .AND. (vbo .LT. 0.0)) THEN 



Tla 


= 


. FALSE . 


T2a 


= 


. TRUE . 


T3a 


= 


.TRUE. 


T4a 


= 


. FALSE . 


T5a 


= 


.FALSE. 


T6a 


= 


.FALSE. 



ELSEIF ((vbo .GT. 0.0) .AND. (vao .GT. 0.0)) THEN 
Tla = .FALSE. 

T2a = .FALSE. 

T3a = .TRUE. 

T4a = .TRUE. 

T5a = .FALSE. 

T6a = . FALSE . 

ELSEIF ((vao . LT . 0.0) .AND. (vco .LT. 0.0)) THEN 
Tla = .FALSE. 

T2a = .FALSE. 

T3a = .FALSE. 

T4a = . TRUE . 

T5a = .TRUE. 

T6a = .FAf.SE 
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ELSEIF { ( vco .GT. 0.0) .AND. (vbo .GT. 0.0)) THEN 
Tla = .FALSE. 

T2a = . FAX.SE . 

T3a = .FALSE. 

T4a = .FALSE. 

T5a = .TRUE. 

T6a = .TRUE. 

ELSEIF ((vbo .LT. 0.0) .AND. (vao .LT. 0.0)) THEN 
Tla = .TRUE. 

T2a = .FALSE. 

T3a = .FALSE. 

T4a = .FALSE. 

T5a = .FALSE. 

T6a = .TRUE. 

ELSEIF ((vao .GT. 0.0) .AND. (vco .GT. 0.0)) THEN 
Tla = .TRUE. 

T2a = .TRUE. 

T3a = .FALSE. 

T4a = .FALSE. 

T5a = .FALSE. 

T6a = .FALSE. 

ENDIF 



END ' " End of derivative" 



DISCRETE intona 

flaga = .FALSE. 

END ! "of discrete inton" 

DISCRETE intonb 

flagb = .FALSE. 

END ! "of discrete inton" 

DISCRETE intone 

flagc = .FALSE. 

END ! "of discrete inton" 

END ! "End of dynamic section" 
END ! "End of program" 
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APPENDIX C. [ACSL PLOTS FOR NBC FAN MODEL WITH 180° INVERTER] 



Machine parameters are those from Table 5.2. 





O 

O 




T »10"3 
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APPENDIX D. [ACSL PLOTS FOR NBC FAN MODEL WITH 120° INVERTER] 



Machine parameters are those from Table 5.1. 
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T * 10~3 
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APPENDIX E. [MATLAB CODE FOR ULTRA CAPACITOR MODEL] 



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%Capt Gabriel Beltran 
%18 Jan 00 

%Ultra Capacitor Model 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% Vendor Data 



figure (1) 

freq = [.01, .0159, .0251, .0398, .0631, .1, .1585, .2512, .3981, . 6310, 1 . O', 1 . 584 9 
2.5119,3. 9811, 6.3096, 10.0,15.8489,25.1189, 39.8107,63.0957. . . 

100.0, 158.4 893, 251.1887,398.1072, 630.9574,1000.0]; 

mag = [23.4411, 14 . 8375, 9 . 6032 , 6. 1331, 3 . 9368, 2 . 5595, 1 . 664 1, 1 . 0782, 0 . 7172 . . 
0.4805,0.3291,0.2325,0.1692,0.1275,0.0991, 0.0787,0.0639,0.0528. . . 
0.0446, 0.0384,0.0337, 0.0300,0.0270,0.0246,0.0226,0.0210] ; 



% Plot 

semilogx (freq, mag) 

title (* Impedance vs. Frequency*) 
xlabel ( * Frequency (Hz) * ) 
ylabel ( * Impedance (Ohms ) * ) 

% Transfer Function Simulation 

figure (2) 
a = 5.350; 

b=input (* input b*) 

num = [b, b*a] ; 

den = [1.0, 0.0]; 

magus = abs (freqs (num, den, freq) ) ; 

% Plot 

semilogx ( freq, magus ) 
title (* Impedance vs. Frequency*) 
xlabel ( * Frequency (Hz) * ) 
ylabel (* Impedance (Ohms) *) 

% Overlay Plots 

figure (3) 

semilogx (freq, mag, * * ’ , freq, magus, * .-* ) 

title (* Impedance vs. Frequency*) 

xlabel ( ’ Frequency (Hz) * ) 

ylabel { * Impedance (Ohms ) ’ ) 

gtext ( ’Actual = ***') 

gtext ( ’Model = .-.-.-*) 
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APPENDIX F. [MANUFACTURER DATA SHEET FOR OPTIMA 800U] 



mm 

^fattery Model: 800U 
'it number: 8004-003 
minal Voltage: 12 volts 
kmi-374-2243 

\ escription : High power ; sealed , lead acid, engine starting battery. 

Ip 

sical Characteristics: 




High purify lead-tin alloy. Wound cell configuration using utilizing proprietary 
SPIRALCELL technology. 



m 



Sulfuric acid, H2S04. Specific gravity. 1.286 



ase: 



Polypropylene 



, 

M 



Case: Dark Gray 
Cover: "Optima" Red 



'mp Size: 



BCI: 34/78 



liillPr • 


Standard 


Metric 


Blfe 


1 C r 


254 mm 


IBS’ 


6.8" 


172.2 mm 




7.8" 


198. 1 mm(heigh( at the top of the terminals) 


WfeigM: 


38.81b. 


17.6 kg 



erminal Configuration: SAE /BCI automotive and GIW style side terminal (318 ■ 16 UNC - 2B, 



i^dednut). 

rfbrmance Data: 



pen Circuit Voltage (full charged): 



jemal resistance (fully charged): 



Reserve Capacity. 



12.8 volts 



0032 ohms 



56 Ah (C/20) 50 Ah (C/2) 



BCI: 120 minutes 

(25amp discharge, 80°F (26.7°C), to 10.5 volts cut-off) 



119 



rer; 



urn: 



M '.(BCI37F): 




mmr(C/2) 




800 amps 



1000 amps 



10.2 kW 

569.8 mg 
1159WIL 



600 Wh 
33.5 Wh/kg 
68.18 Wh/L 



lUSe life - deep cycle (BCI): 


~50 cycles 


mM$m Cycle {BCD- 


10,000 -12,000 cycles 



commended Charging: 

Mowing charging methods are recommended to ensure a long battery fife: 
ays use a voltage regulated charger, with voltage limits set as described below.) 

?od# 800U 

ie batteries are designed for engine starting applications. They are not recommended or 
Tainted for use in deep cycle applications. 

Recommended charging information: 



13.8 to 15.0 volts no amperage limit 



13.8 to 15.0 volts, 10 amps maximum, 8 hours maximum 



K 

ost Recharge 



Maximum voltage: 15.6 volts (regulated). 
Maximum current: 60 amps. 

Maximum recharge time: 2 hours. 

All limits must be strictly adhered to 



13.2 to 13.8 volts, 1 amp maximum current, time indefinite (at lower voltage levels). 
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toys wear safety glasses when working with batteries. 

i rays use a voltage regulated battery charger with limits set to the above ratings. Overcharging 
i cause the safety valves to open and battery gasses to escape, causing premature failure. These 



Mrs 



iy battery that becomes very hot while charging should be disconnected immediately. 



sir 






( fully charging a battery can result in poor performance and a reduction in capacity. 

■Site 



ipping and Transportation information : 



Ew*' , 



r Sima batteries can be shipped by AIR. The battery is non-spillable and is tested according to 
{ CAf^jMnical Instructions DOC. 9284-AN/905 to meet the requirements of Packing Instructions 
fig 80S and is classified as non-regulated by IATA Special Provision A48 andA-67 for UN2800. 



RK3£ 



rertnina/s must be protected from short circuit. 



j during Location : 



vma Batteries, Inc. 
~'JEast 22nd A venue 
CO 80011 

fhited States of America 
ie: 303-340-7400 
103-340-7474 



Jg Battery Council International 
^European Norm 



litharge versus temperature (800 amps): 
"Marge versus temperature (400 amps): 
charge: 

Sneering Drawing: 

pip 

Ana Batteries, Inc. 
futt Specifications: Model 800U 
7/0M7 
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